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Abstract 

A systematic study of electronic structure and band gap states is conducted to 
analyze the mono doping and charge compensated codoping of rutile Ti02 with Rh 
and Nb, using the DFT+J7 approach. Doping of rutile Ti02 with Rh atom induces 
hybridized O 2p and Rh 4d band gap states leading to a red shift of the optical ab- 
sorption edge, consistent with previous experimental studies. Since Rh mono-doping 
may induce recombination centers, charge-compensated codoping with Rh and Nb is 
also explored. This codoping induces an electron transfer from Nb induced states to 
Rh 4d states which suppresses the formation of Rh"*""", thereby leading to a reduction 
in recombination centers and to the formation of more stable Rh^^. A combination 
of band gap reduction by 0.5 eV and the elimination of band gap states that account 
for recombination centers makes (Rh,Nb)- codoped TiOi a more efficient and stable 
photocatalyst. 



1 Introduction 

The development of a safe, abundant and inexpensive photocatalytic material is an impor- 
tant consideration in solving our imminent energy and environmental crisis [1]. Titanium 
dioxide (Ti02) is one such promising material that meets the requirements of a good 
photocatalyst and has therefore attracted wide attention in recent years for use into solar 
hydrogen production [Q,!!], water decontamination [0], and solar cells [H]. However, TiOa 
is a wide band gap semiconductor (3.0 eV for rutile and 3.2 eV for anatase). This restricts 
sunlight absorbance in the UV range, accounting for < 5% of the full solar spectrum and 
rendering pure Ti02 inefficient for solar energy conversion. Suitable reduction of the band 
gap, either by doping or by surface engineering, can lead to significant improvements in 
the photocatalytic activity of Ti02 [@,l3,[8l,@]- Doping leads to defect states in the band 
gap [jlOjl which allow sub-band-gap illumination by electronic transitions from VB to de- 
fect states or from defect states to CB thereby improving visible light activity. However, 
these intraband states may also serve as recombination centers leading to a significant 



drop in the overall photocatalytic efficiency nll\\ 



Ti02 doped with rare earth metals possess significantly improved photo activity owing 
to the enhanced electron density that is imparted to the titania surface by the dopant oxide 



UM- Despite their potential significance, investigations on doping with rare earth metals 
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have so far been limited. In particular, experimental reports of rhodium (Rh) doping 
suggest that photocatalytic activity depends on the dop ant concentration, its distribution, 
visible light intensity, and the preparation method [|lll.ll3ll. Choi et al. [Illll reported that 
doping Ti02 with Rh increases the photo reactivity markedly, due to the charge transfer 
transitions between Rh and Ti02 CB or VB; or a d-d transition in the crystal field. Kitano 
et al. ids'] showed that Rh modified Ti02 samples exhibit higher level of activity than 
Cu^"^ modified and the nitrogen-doped samples. On the other hand, Niishiro et al. [14] 
reported that Ti02 doped with only Rh shows only negligible photocatalytic activity as 
it creates Rh"*"^ ions. They suggested that in such cases, charge compensated codoping 
with pentavalent ions such as Sb^^, Nb^^, Ta^"^ becomes necessary for suppressing the 
formation of Rh"^^ ions and improving photocatalytic activity. 

Motivated by above experimental findings and apparent discrepancy in the role of 
Rh-dopant, we performed a detailed investigation of the band structure and electronic 
properties of Rh-doped, Nb-doped and (Rh,Nb)- codoped rutile. Codoping was investi- 
gated in order to understand the potential benefit of charge compensation for eliminating 
or suppressing the recombination centers that may be present in mono-doped rutile struc- 
tures. First principles based theoretical investigations were conducted so as to clearly 
delineate the origin of band gap defect states, the local structural environment around the 
dopant site and the localization of excess charge due to defects. Furthermore, in order 
to represent strong d-electrons correlation in titania, the present study utilized the density 
functional theory with Hubbard energy correction, i.e. the DFT+U methodology. An in- 
depth comparison of our theoretical findings with above mentioned experimental data is 
performed so as to ascertain the effect of rare earth metal doping of rutile in the hope of 
developing an efficient visible light photocatalyst. 



2 Computational details 

Rutile has a tetragonal crystal structure, with a space group of p42/mnm containing two 
Ti and four O atoms. A 2x2x3 supercell consisting of 72 atoms was considered in the 
present calculations. The Quantum-Espresso flS*] software package within the plane- 
wave pseudopotential approach was utilized throughout. The interaction between the 
valence electrons and the ionic core was described by the generalized gradient approx- 



imation (GGA) in the Perdew-Burke-Emzerhof (PBE) formulation 111611 with Vanderbilt 
ultrasoft pseudo-potentials [IitIi . The kinetic energy cutoffs of 544 eV (=40 Ry) and 5440 
eV were used for wave functions and charge density, respectively. Brillouin zone inte- 
grations were performed using a Monkhorst-Pack [18] grid of 4 x 4 x 6 k-po'mts. All 
calculations were spin polarized. The structures were relaxed using a conjugate gradi- 
ent minimization algorithm until the magnitude of residual Hellman-Feynman force on 
each atom was less than 0.025 eV/A. The bulk lattice parameters for zero pressure rutile 
structure obtained by variable cell relaxation using DFT, a=b=4.615 A and c=2.961 A, 



were found to be in good agreement with experiments [|l9n and other theoretical calcula- 
tions H]. 

Appreciable underestimation of band gap and delocalization of d and f electrons are 
well known limitations of DFT. Therefore, we used DFT-l-i7 formalism in these calcu- 
lations. Our DFT calculations suggested a band gap value of 1.89 eV, consistent with 
other theoretical calculations [|2li1 . On the other hand, DFT+U approach with U applied 



2 



to Ti 3d electrons, varying from 1 to 11 eV, revealed a band gap of 2.85 eV at U = 8 
eV, far closer to the experimental value. Nevertheless, choosing U solely on the basis of 
band gap may lead to incorrect conclusions regarding the effect of dopants on the elec- 
tronic properties ll22ll . It must rather be chosen depending on the property of interest [|23ll . 
which, in the current study, is the photocatalytic behavior of Ti02 that in-tum depends 
upon the position of band gap states and their effect on the electronic structure. Spectro- 
scopic experiments for bulk reduced rutile Ti02 crystals suggest a broad gap state lying 
about 0.9 eV below the CB edge JmIiSD. Following the approach of MoUer et al. [26], it 
was found that U= 4.2 eV yields correct position of band gap states due to Ti interstitial 
and O vacancy defects, while opening up the band gap to 2.2 eV for rutile. This value 
of U is consistent with theoretical investigations by Morgan et al. ^T^, who calculated 
it by fitting the peak positions for surface oxygen vacancies to experimental XPS data. 
Furthermore, the resulting electronic projected density of states (PDOS) plot for Ti inter- 
stitial and O vacancy shows that these defect states are composed of mainly titanium d 
states, in good agreement with experiments II28I1 . 



3 Results and discussion 



Experimentally obtained X-ray diffraction peaks of Rh-doped TiOi have revealed that 



Rh is substituted for Ti sites 



This is expected because Rh^^ ions (0.805 A) are 
bigger than Ti'*^ ions (0.61 A) [|29(]. As radius of Nb (0.64 A) is also comparable to Ti 
atom, it should also be substituted to Ti atoms. The concentration of Rh and Nb ions in 
present study is ~1.4% commensurate with the concentration of dopants considered in 
experiments nl4\ . 



3.1 Structural distortion 

Doping leads to local structural disorder that modifies the electronic environment and 
therefore affects photocatalytic properties. Here, structural distortion was investigated by 
analyzing the local coordination environment, nature of chemical bonding, and dopant-0- 
dopant distances. Undoped rutile has a symmetric structure where each Ti atom is bonded 
to four nearest and two second-nearest oxygen neighbors with bond lengths of 1 .95 A and 
1.99 A, respectively. Introduction of an Rh dopant into rutile leads to the elongation 
of all six Rh-0 distances without breaking symmetry. However, introduction of the Nb 
dopant leads to elongation of four Nb-0 distance and shortening of the other two Nb-0 
bonds, thereby affecting the symmetry around the dopants to an appreciable extent. A 
comparison of ionic radii of dopants with respect to Ti atoms suggests that a dopant with 
smaller ionic radius leads to larger structural distortion. The local structure around the 
Rh and Nb in codoped Ti02 is strongly distorted as compared to Rh and Nb monodoped 
Ti02. For (Rh,Nb)- codoped structure the 0-Rh bonds lengths are 2.038, 2.063, and 2.029 
A and 0-Nb bond lengths are 2.005, 1.997 and 1.925 A. Thus, in (Rh,Nb)- codoping four 
Nb-0 distances get elongated, two Nb-0 distances get shortened, and all 6 Rh-0 distances 
get elongated. Overall, Nb^^ attempts to shorten two cation-0 distances and leaves most 
of the Rh-0 distances notably elongated in order to obtain a coordination environment 
favorable for its ionic radius. This results in a comparatively more distorted structure than 
mono doped Ti02 structures. 
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3.2 Formation energy 



The formation energy is an important criterion to evaluate the relative difficulty for the 
incorporation of dopants into the host lattice and therefore the relative stability of the 
doped systems. The formation energies of Rh-doping, Nb-doping and (Rh,Nb)- codoping 
can be written as 

EfiRh) = E{Ti24~niRhni04g) - (£(r?24048) + RiyUj - HoUo), (1) 

Ef(Nb) = E{Ti2A-„^Nbn^04^) - {E{Ti2A0^^) + n2jU2 - hqHq), (2) 

Ef(RhNb) = E(Ti24^ni-n2RhrnNbn204s) - (,E(Ti2404s) + TliHi + n2H2 - noflo), (3) 
where E{Ti24-n^Rhn,04^), E{Ti24-n2Nbn204&), E{Ti24-n,-n2RhniNbn204%) and E{Ti2404%) 

are the total energies of the Rh-doped Ti02, Nb-doped Ti02, (Rh,Nb)- codoped Ti02 
and pure Ti02, respectively, hq, jxi and H2 are the chemical potentials of the Ti, Rh and 
Nb elements, respectively; and hq, n\ and n2 denote the numbers of the host (Ti) atoms 
substituted, Rh-dopant atoms and Nb-dopant atoms, respectively. In the present study, 
hep bulk metal Ti, fee bulk metal Rh and bcc bulk metal Nb were used to determine 
the chemical potentials. Since the focus here is on the relative stabilities, the choice of 
chemical potential should not influence the conclusions about stability of the doped Ti02. 
The calculated formation energies for Rh-doped and Nb-doped rutile were -1-0.417 eV 
and -0.0689 eV, respectively. Postive and negative values of formation energy for Rh- and 
Nb-doped Ti02 respectively, suggest that Nb-doped Ti02 is more stable as compared to 
Rh-doped Ti02. 

Calculation of formation energy is also useful in determining the most stable configu- 
ration of the codoped structure. For this purpose, six nearest-neighbor codoping systems, 
labeled 1-6 in Fig. [U were constructed. The calculated formation energies for Rh-Nb 
distances of 2.95, 3.68, 4.62, 5.56, 5.63 and 7.23 A were respectively equal to 0.232, 
0.241, 0.240, 0.239, 0.240 and 0.233 eV. A relatively smaU range of formation ener- 
gies of structures with difl^erent Rh-Nb distances indicates high stability of the codoped 
Ti02 structures. For subsequent study of electronic properties, the most stable (Rh,Nb)- 
codoped structure with dopant-dopant distance of 2.95 A and Ef=0.232 eV was cho- 
sen. The chosen Rh and Nb atoms are highlighted in blue and green respectively, in 
Fig. [TJ (Rh,Nb)- codoped Ti02 has a smaller formation energy than that of Rh-doped 
Ti02. Therefore, choosing Nb as codopant not only favors charge compensation but may 
also lead to spontaneous formation of (Rh,Nb)- codoped system. This finding is consis- 
tent with experimental observations [14]. On the whole, structural distortion and forma- 
tion energy results suggest that stability of defect formation depends on the ionic radius of 
the dopant, in agreement with Nolan 113011 . The strong distortion due to small ionic radius 
of dopant Nb makes (Rh,Nb)- codoped Ti02 more stable than Rh-doped Ti02. 

3.3 Electronic properties 

In order to understand the photocatalytic properties of doped systems, the spin polarized 
density of states (DOS) and PDOS were calculated for Rh-doped, Nb-doped and (Rh,Nb)- 
codoped rutile Ti02. The corresponding spin density plots are also shown. A conventional 
Gaussian smearing of 0.1 eV was utilized. 
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3.3.1 Rh-doped rutile TiOa 



The total DOS with PDOS for the d electrons of Rh and Ti atoms and p electrons of 
O atoms of the Rh-doped rutile are summarized in Fig. Oa). Rh doping increases the 
asymmetry between the neighbouring atoms and introduces an unpaired d electron from 
Rh. As a result an intermediate energy band of 0.3 eV width is found approximately 1.0 
eV below the CB minima, along with some states at CB and VB tails. This leads to an 
overall reduction in the band gap by 0.5 eV. To understand the nature of defect states 
further, consider the electronic configurations ofRh: 4d^5s^ and Ti: 3d^4s^. Substitution 
of Ti with Rh leads to a formal +4 state (Rh'^^), i.e. 4d^5s" configuration. Since the stable 
oxidation state of Rh is +3, an electron get transferred from neighbouring O atoms to 
Rh'*^ leading to O^^ defects. The band gap states are mainly due to hybridization of Rh 
4d and O 2p states. The spin density plot (Fig. [2tb)) also shows that the spin densities 
are strongly localized on the dopant atom and O'" sites in its neighborhood. To further 
clarify the origin of band gap states, we present in Fig. [2c) the five-fold degenerate Rh-4d 
states. The 2-fold degenerate e^ states split into dji2_y2 and d^2 states; whereas the 3-fold 
degenerate t2g states split into doubly degenerate states (d^z, dv,) and a single state (dx,.). 
It is evident from Fig. [2tc) that only dxi^yi state of Rh-4d orbital contributes to band gap 
state. Rh doping of rutile has also been investigated by Tan et al. yjj]. Their calculations 
suggested band gap states very near the VB maxima, whereas in the present case defect 
states are almost in the middle of band gap. Nonetheless, in both studies, the contribution 
to band gap states comes from the same d,-2_^,2 state of Rh-4d orbital. The difference in 
position of band gap states can be attributed to different choice of U parameter. 

The formation of a VB by orbitals not associated with O 2p but with other elements is 
indispensable in designing visible-light-driven photocatalysts as it moves top of the VB 
at a more positive potential than the oxidation potential of H2O to O2 without affecting 
the CB level and charge migration i32ll . Thus an acceptor state at the VB maxima not 
associated with O 2p but with the Rh"^^ ions can increase photo activity in Rh-doped TiOi 
significantly. Additional states present in the band gap may also contribute to increased 
photo activity due to electron excitations under visible light irradiation from top of the 
VB to the band gap states and then from band gap states to the bottom of the CB. Choi et 
al. [flUl also suggests that Rh-doped Ti02 shows red shift due to charge transfer between 
Rh and Ti02 CB or VB; or a d-d transition in the crystal field. These mid-gap states 
are not without a downside, however, as they can act as recombination centers, conse- 
quently leading to a significant reduction in the photocatalytic efficiency. As per Niishiro 
et al. [14], this recombination is chiefly due to the formation of Rh'^^. To suppress Rh'*^, 
codoping of Rh-doped rutile with high-valent Nb^^ could be useful in limiting forma- 
tion of these recombination centers, and is therefore investigated further in subsequent 
discussion. 



3.3.2 Nb-doped rutile Ti02 

Before proceeding to charge compensated doping, let us describe observations on Nb- 
monodoped rutile. From the total DOS, PDOS and accompanying spin density plots 
shown in Fig. [31 it can be deduced that the substitution of Ti"*^ atom with Nb^^ in Ti02 
results in a narrow peak of 0.3 eV width in the band gap approximately 1.3 eV below 
the CB edge. It is noted here that our independent calculations using DFT alone did not 
show a band gap state with excess charge occupying the bottom of the CB. This implies 
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the necessity and significance of the DFT+t/ approach in reproducing electronic states 
and band structure commensurate with the experimental observations in an accurate and 
physically consistent manner. The PDOS plot (Fig.[3ta)) also suggests that the defect state 
is mainly due to contribution from Ti atoms. Experimental data have also shown a defect 
state in the band gap well separated from the CB [[330, with a height proportional to the Nb 
concentration. This observations on Nb-doped Ti02 with DFT+U agrees with previous 
theoretical investigations by Morgan et al. 1340 ■ Considering the electronic configuration 
of Nb: 4p^4d^5s^, substitution of Ti by Nb should resuh into Nb"^"" state, i.e. 4p^4d'5s" 
configuration. To obtain a stable oxidation state of +5 the extra electron from Nb'*^ is 
preferentially transferred to Ti'*^ and therefore Ti^^ state is formed. This is confirmed 
by spin density plot (Fig Ob)) which shows that the spin density strongly localized on a 
single Ti site in a nearest-neighbour position relative to the Nb dopant. 

3.3.3 (Rh,Nb)- codoped rutile TiOz 

The presence of Rh states at VB maxima and intermediate band with considerable band- 
width in Rh-doped rutile indicates great potential of using the Rh-doped material for vis- 
ible light photocatalytic applications. This intermediate band offers a stepping stone for 
the absorption of low energy photons via the excitation of electrons at the VB maximum 
to the intermediate bands, from where they can be excited again above the CB for effective 
photoactivity of Ti02. Nonetheless, these band gap states can also act as recombination 
centers, see e.g. experimental observations by Niishiro et al. [14]. This charge recombi- 
nation renders Rh-doped Ti02 as ineffective photocatalyst. These recombination centers 
need to be eliminated by suppressing Rh'*^ ions. This can be achieved by codoping of 
TiOi by a combination of p-type and n-type dopants, as also suggested by [fl4[ BSll. In 
this context, further codoping of Rh-doped rutile with Nb was investigated with the goal 
of eliminating of Rh"^^ ions through charge compensation. For studying codoped system, 
we constructed 6 nearest neighbour codoped structures. These structure were relaxed and 
the structure having Rh-Nb distance of 2.95 A showed lowest configuration energy and is 
therefore considered for further analysis. 

Total spin polarized DOS and PDOS for the d electrons of Rh, Nb and Ti atoms and 
p electrons of O atoms of the (Rh,Nb)- codoped rutile are shown in Fig. Ufa). Based on 
the analyses of Rh and Nb monodoping described earlier in the paper, Rh dopant requires 
an extra electron to achieve the stable -1-3 oxidation state whereas Nb dopant introduces 
an extra electron into the system. Thus, when rutile is codoped with Rh and Nb the ex- 
tra electron on Ti^^ due to Nb transfers to Rh dopant such that the Rh species attains a 
-1-3 state (4d^5s') while the Nb species attains a +5 (4p^4d°5s") state, with both result- 
ing states being more stable. This charge compensation of (Rh,Nb)- codoped system can 
also be inferred in the spin density plot (Fig Sfb)) by the absence of excess spin density 
states. This codoping also reduces the band gap by 0.5 eV relative to undoped Ti02, due 
to the new VB consisting of hybridized Rh and O 2p orbitals which results in decreased 
band gap energy without affecting the CB level. The Rh d^2_Y2 orbital, which was earlier 
responsible for band gap states in Rh-doped rutile, now shifts to a fully filled state at va- 
lence band minimum, as shown in Fig.H^c). Furthermore, the codoped structure does not 
require formation of O vacancies to provide extra electron to Rh. This, in turn leads to 
a smaller degree of defectiveness and elimination of recombination centers. Our analy- 
sis thus establishes the beneficial aspects of charge compensated codoping and seems to 
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explain findings of Niishiro et al. [14] for (Rh,Nb)- codoped Ti02. 

3.4 Charge density 

To substantiate the electronic properties further and to study the variation of chemical 
bonding induced by Rh and Nb codoping, the electron charge density distributions for 
undoped and charge compensated systems were also analyzed. As depicted in Fig. [51 the 
contour plots of the charge density on (1 10) and (0 1 0) represent the effect of dopants 
on nearest neighbouring O and Ti atoms, respectively. The charge density around Rh and 
Nb atoms is nearly spherically distributed with a slight deformation toward their nearest- 
neighboring atoms. Fig.[5te) shows that there are clear covalent bridges between dopants 
and O atoms, which strengthens the judgement of certain covalent character of (Rh,Nb)- 
codoped rutile. Furthermore, it can be seen that the overlap of Nb 4d orbital with O 2p is 
greater than that of Rh 4d orbital with O 2p, which indicates that the degree of covalency 
between Nb and O is stronger than that between Rh and O. A weaker interaction between 
Rh and O atoms as compared to that between Ti and O atoms makes the Rh-0 bond longer 
and leads to smaller binding energy. Consequently the O 2p states, which are dominating 
in valence band edge, move toward high energy levels and shift the VB maximum to 
increased energy, thereby reducing the band gap. From (0 1 0) plane, on comparing with 
undoped rutile (Fig. [2d)), it can be seen that there is a considerable charge accumulation 
in the bonding regions of Rh and Nb dopants (Fig. [3f)). This confirms the sharing (or 
transfer) of electron between the dopants. The accumulation of charge in the region close 
to the dopants also suggests the formation of a metallic contact. It should be noted that 
the electron charge density plotted on plane next to doped elements is almost completely 
unaffected by the presence of the doping elements. 

4 Conclusions 

In this work, ab-initio calculations were conducted to understand conflicting experimen- 
tal results for photocatalytic behavior of Rh-doped Ti02 and to explore a more efficient 
photocatalysts by codoping it with Nb. DFT+U methodology was utilized to remove ex- 
cessive delocalization of Ti d-orbitals present in standard DFT calculations. Doping rutile 
with Rh atom induces states at valence band minima and an intermediate band that leads 
to a red shift in optical absorption edge; also observed in previous experimental inves- 
tigations. But since these states may also act as recombination centers due to presence 
of Rh"*^ ions, codoping with Nb was investigated. Following conclusions can be derived 
from the present study: 

1 . The formation energy of Rh-doped system reduces appreciably on codoping it with 
Nb, suggesting a spontaneous formation of (Rh,Nb)- codoped rutile Ti02 at room 
temperature. 

2. Codoping affects the symmetry, leading to a weaker interaction and longer bond 
between Rh and O as compared to Nb and O. This structural distortion reduces the 
binding energy and the O 2p states tend to move from VB edge towards high energy 
levels, reducing the overall band gap. The hybridized O 2p and Rh 4d states at the 
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top of VB in codoped rutile are thought to be responsible for the red shift of optical 
absorption edge, confirming previous experimental studies [|l4 llU ll3n. 



3. Rh dopant requires an extra electron to achieve the stable +3 oxidation state whereas 
Nb dopant introduces an extra electron into the system. Thus, when the system is 
codoped with Rh and Nb, both dopants attain their stable oxidation states. It means 
that the creation of oxygen vacancy is not required in a codoped structure. This 
leads to a significant reduction in the level of defectiveness and thus an enhanced 
photocatalytic activity. 

4. Unlike Rh-doped rutile that has an isolated band gap state which may act as a re- 
combination center, charge compensated (Rh,Nb)-codoping does not show such 
bandgap states which commensurate with experimental findings of Niishiro et 
al. 1 14]. This elimination of recombination centers in combination with a band 
gap reduction by 0.5 eV makes codoped rutile a better photocatalyst. 

Overall, the present study reinforces the usefulness of charge compensated codoping for 
reducing the band gap in photocatalyts, without creation of unwanted recombination cen- 
ters. In the near future, we plan to apply this technique further to amorphous Ti02 in the 
hope of developing more efficient and cost-effective visible light photocatalysts. 
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Figure 1 : The (Rh,Nb)- codoped 2x2x3 supercell of rutile. Ti, O, Rh and Nb atoms are 
highlighted in gray, red, blue and green colors, respectively. The six nearest neighbouring 
positions considered for Nb atom are shown along with dopant-dopant distances. 
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Figure 2: The density of states and projected density of states on the p and d orbitals (a), 
spin density (b) and projected density of states for the Rh-4d states (c) in Rh-doped rutile 
model. The zero energy value is set at the top of VB and Fermi energy is represented by 
the vertical dashed line in (a) and (c). 
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Figure 3: The density of states and projected density of states on the p and d orbitals 
(a) and spin density plot (b) for one Nb atom doped in 72-atoms rutile model. The zero 
energy value is set at the top of VB and Fermi energy is represented by the vertical dashed 
line in (a). 
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Figure 4: The density of states and projected density of states on the p and d orbitals 
(a), spin density (b) and projected density of states for the Rh-4d states (c) in (Rh,Nb)- 
codoped rutile model. The zero energy value is set at the top of VB and Fermi energy is 
represented by the vertical dashed line in (a) and (c). 
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(a) (b) 




Figure 5: Two dimensional 1 10 (a) and 010 (b) planes indicated in black for which charge 
density contour plots (electron bohr"^, (lbohr= 0.529 A) units) are depicted. The charge 
density plot for undoped rutile in 1 10 plane (a), undoped rutile in 010 plane (b), (Rh,Nb)- 
codoped rutile in 110 plane (c) and (Rh,Nb)- codoped rutile in 010 plane (d). The red 
regions denote the area with nearly zero electron density and the pink regions denote the 
highest electron density. 
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